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Recent realizations of single-atom trapping and tracking in cavity QED open the door for feedback 
schemes which actively stabilize the motion of a single atom in real time. We present feedback 
algorithms for cooling the radial component of motion for a single atom trapped by strong coupling 
to single-photon fields in an optical cavity. Performance of various algorithms is studied through 
simulations of single-atom trajectories, with full dynamical and measurement noise included. Closed 
loop feedback algorithms compare favorably to open-loop "switching" analogs, demonstrating the 
importance of applying actual position information in real time. The high optical information rate 
in current experiments enables real-time tracking that approaches the standard quantum limit for 
broadband position measurements, suggesting that realistic active feedback schemes may reach a 
regime where measurement backaction appreciably alters the motional dynamics. 



I. INTRODUCTION 

Recent experiments |B| have demonstrated 

the ability to trap 0, 1111^13^ and localize @ a single 
atom in a high-finesse optical cavity by way of optical 
forces. Moreover, by detecting the light transmitted by 
the cavity, the atom's motion within the cavity mode 
can be monitored in real time with high signal to noise 
throughout its trapped lifetime 7]. These achievements in 
trapping and localization open up exciting possibilities 
for quantum logic and quantum state preparation in the 
context of cavity QED|1 U [H El El El Beyond the 
realization of trapping, the high signal-to-noise for con- 
tinuous, real-time position measurement is itself one of 
the most notable features of these strongly coupled cavity 
QED systems. Such detailed real-time position informa- 
tion immediately suggests the idea of active feedback to 
dynamically cool the motion of a single trapped atom. 
By investigating this system and the feedback schemes 
available in it, basic questions of quantum state estima- 
tion and optimal control can be explored for continuous 
measurement of a dynamical variable - in this case the 
position of a single atom. 

Crucial to the realization of trapping and sensing in 
cavity QED is strong coupling, a condition in which the 
coherent coupling between atom and cavity field dom- 
inates dissipative rates in the system. For a two-state 
atom optimally coupled to the cavity mode, the dipole- 
field coupling is given by the Jaynes-Cummings interac- 
tion Hamiltonian [l4| 

H lnt = hg(a + a + er-a 1 ), (1) 

where a± are dipole raising and lowering operators, 
(a,a^) are field annihilation and creation operators for 
the cavity mode, and g is one half of the single-photon 
Rabi frequency. This interaction gives rise to the well- 
known Jaynes-Cummings ladder of eigenstates for the 
coupled atom-cavity system, and correspondingly to the 
vacuum Rabi splitting for the system's resonant frequen- 
cies ^3 • Dissipation, on the other hand, is characterized 
by the cavity decay rate n and the atomic spontaneous 
emission rate 7. Strong coupling occurs for g 3> (k, 7). 



We can define a further condition of strong coupling for 
the external atomic degrees of freedom; this occurs when 
the coherent coupling also dominates the atomic kinetic 
energy, as first achieved in Ref. ^3 . 

Under these conditions, interaction with a single- 
photon cavity field exerts a strong mechanical effect on a 
single atom, allowing trapping of the atom when the sys- 
tem is driven to strong-field seeking states in the Jaynes- 
Cummings ladder. Furthermore, strong coupling assures 
that the intracavity light field and thus the cavity output 
field (transmitted light) are influenced by an atom; thus, 
as an atom moves between more and less strongly coupled 
positions within the cavity, the transmitted light pro- 
vides a real-time measurement of atomic position. This 
sensing enables the one-time triggering employed in Refs. 
IE tl t° switch on a trapping potential when an 
atom is present near the center of the cavity. The ongo- 
ing stream of position information should, however, be 
useful for continued active feedback based on the atomic 
position. Rather than simply triggering a potential to 
turn on, it should be possible to modulate the poten- 
tial depth to dynamically cool the atomic motion, in a 
method analogous to the principles of stochastic cooling 
but for a single atom. Initial steps in this direction are 
represented in the work of Ref. [TtJ ■ Quantum feedback 
for other atomic state variables is also an active area of 
research, including the recent experimental demonstra- 
tion of quantum feedback to control the ensemble spin of 
a collection of cold atoms [Is| . 

In this paper, we address the question of how to im- 
plement atomic position feedback in experimentally real- 
istic situations, where several constraints apply: (1) the 
system is inherently nonlinear and largely nonanalytic, 
with relationships between many quantities of interest 
determined by steady-state solutions to the master equa- 
tion for the atom-cavity system, (2) dynamical noise is 
significant and changes in tandem with the driving field 
and trapping potential, and (3) measurement noise, aris- 
ing largely from the fundamental quantum noise (shot 
noise) of detection, imposes necessary delays in the es- 
timation of dynamical variables and the implementation 
of feedback. Broadband measurement near the standard 



2 



quantum limit (see Ref. 19]) has been demonstrated in 
this system through measurement of cavity transmission 
amplitude pj and by simultaneous measurement of trans- 
mitted amplitude and phase [2(}. Use of these measure- 
ments in feedback control of some aspect of the atomic 
motion should bring us closer to regimes where measure- 
ment backaction has a significant effect, so that different 
detection methods may exhibit different control limits 
based on these effects as well as on more conventional 
signal-to-noise considerations 0, • 

In Section II we review the experimental system on 
which our work is based. The feedback strategy we con- 
sider is presented and motivated in Section III. Section 
IV presents feedback results; these are considered both in 
the realistic experimental context and in several idealized 
systems in order to illustrate the algorithm's effect on the 
individual components of atomic motion. In Section V 
we again take up the topic of experimentally measurable 
feedback results, and in Section VI we discuss limits and 
possible extensions of the algorithms discussed in the pa- 
per. 



II. EXPERIMENTAL STATUS AND 
SENSITIVITY FOR ATOMIC POSITION 
MEASUREMENT 

In our feedback calculations, we have considered the 
situation of Ref. , where a single atom is trapped via 
its interaction with a near-resonant cavity light field at 
the level of a single photon. In that experiment, the small 
saturation photon number and critical atom number p3 | 
- the smallest in experiments to date - facilitated not 
only the trapping mechanism but, more importantly for 
these considerations, the high signal-to-noise observation 
of atomic motion within the cavity field. 

Figure n shows a schematic diagram of the experimen- 
tal procedure, in which Cesium atoms (atomic resonance 
frequency uj a ) are collected in a magneto-optical trap 
(MOT) directly above the cavity mirrors, cooled to tem- 
peratures of ~ 10/xK, and dropped through the cavity. 
The geometry of the mirror substrates cuts off most of 
the atomic flux so that one atom at a time transits the 
cavity mode of length I and Gaussian waist wq. The 
cavity resonance oj c is tuned near but slightly below the 
atomic resonance frequency so that uj c — uj a < 0. The cav- 
ity is continuously driven by a probe laser at frequency 
top, and the transmission of this beam through the cavity 
is monitored via balanced heterodyne detection. For a 
probe red-detuned from both atom and cavity (uj p near 
the lower vacuum Rabi sideband 15]), transmission is low 
for the empty cavity and is highest when an atom is 
in the regions of strongest coupling. Thus the moni- 
tored photocurrent carries real-time information about 
the atomic position, with high signal-to-noise even for 
probe strengths corresponding to < 1 intracavity pho- 
tons. Saturation of the atom-cavity response|l6| sets in 
for larger field strength, so that the most sensitive track- 



ing is realized at ~ 1 intracavity photon. 




FIG. 1: Cesium atoms are dropped from a MOT and a small 
fraction of these fall one by one through the mode of an optical 
cavity. The cavity mode is a cylindrically symmetric Gaus- 
sian in the transverse (p, 9) plane and has a standing- wave 
structure in the axial (x) direction. Transmission of a probe 
beam through the cavity is monitored to sense an atom's mo- 
tion through the cavity mode; a detected rise in transmission 
triggers a switch to higher probe intensity, significantly pop- 
ulating strong-field-seeking states of the atom-cavity system 
and thus trapping the atom. The figure shows one falling 
atom which misses the cavity mode altogether and a second 
atom which enters the cavity, is trapped in the mode, and 
eventually escapes. 

The atom-cavity evolution in thi s sy stem is described 
by a master equation (see, e.g., [2J, I251 |26j ) for the joint 
atom-cavity density operator p. We consider a driving 
(and probing) field e of frequency u) p , a cavity resonant 
at oj c = uj p + A cp , and an atomic resonance frequency 
uj a = uj p + A ap . In the electric dipole and rotating- wave 
approximations, and in the interaction picture with re- 
spect to the probe frequency, the master equation can be 
written 

p = -~[H ,p] + 7(2(Tp<7 t - a^ap - pa^a) 

-\-K,(2apa) — a)ap — per a), (2) 
Hq = HAcpC^a + hA ap a^a + hg{r)[aa^ + a) a] 

+fie(a + a t ). (3) 

Here g(r) is the coupling strength which takes into 
account the atomic position r within the cavity 
mode. For a Fabry-Perot cavity supporting a stand- 
ing wave mode with Gaussian transverse profile, g(f) = 
gocos(2irx / A)exp[— (y 2 + z 2 )/w 2 ]. The cylindrical sym- 
metry of the field suggests the use of cylindrical co- 
ordinates (p,6, x), where p = \J y 2 + z 2 and 9 = 
tan~ 1 (—z/y) (see Figure QJ. Thus we write g(r) = 
gocos(2irx/\)exp[—p 2 /w 2 ,}. In the fully quantum treat- 
ment, the atomic position f is itself an operator; in ex- 
periments to date a quasi-classical treatment suffices, so 



the atom may be considered a wavepacket with r a clas- 
sical center-of-mass position variable. Similar feedback 
schemes for an atom-cavity system have also been ex- 
plored theoretically in the case of an atom which has 
already been cooled radially and must now be treated 
in a fully quantized manner for cooling of the remaining 
axial motion [23. 

Following the experimental situation of Ref . pj , we con- 
sider an atom-cavity system in which (go, k,j)/2t: = 
(110,14.2,2.6) MHz. The simulation results below re- 
fer to varying cavity field strength but with detunings 
fixed at (w c -w a )/27r = -47 MHz, (w p - uj a )/2n = -125 
MHz. 

Figure OJa) shows the first few levels of the Jaynes- 
Cummings ladder of energy eigenvalues, obtained by di- 
agonalizing the interaction Hamiltonian of Eq. [3J The 
smooth evolution from uncoupled to fully coupled eigen- 
states reflects the dependence of coupling g(r) on the 
atomic position r, and specifically on the atom's distance 
p from the cavity axis. In the presence of dissipation 
and driving, the distribution of populations across the 
first few levels of this ladder is determined by numeri- 
cal steady-state solution of the master equation at each 
position. Figure |2Jb) shows the effective potential for 
the atom-cavity properties considered in this work, with 
driving level in this example fixed at rihi—0.3 photons 
in the empty cavity. The atom-trapping scheme of [l| is 
based on tracking the atomic position and altering the 
driving field strength to place the system in the attrac- 
tive potential of Figure [^b) when the atom is close to 
the cavity axis (p~0). 

Figure EJa) shows a sample experimental trace of 
transmission vs. time for a single atom trapped as in 
Ref. Q . Immediately notable in the transmission record 
are the large, regular oscillations in the heterodyne cur- 
rent, which can be associated with an atom repeatedly 
approaching and receding from the regions of strongest 
atom-cavity coupling. The cavity mode structure pro- 
duces a trapping potential with width ~ A/ 2 = 426nm 
in the axial direction and ~ wq = 14/im in the radial 
direction, giving axial oscillations at ~ 1 — 2MHz while 
radial oscillations occur on the much slower timescale of 
~ 10— 20kHz. By using a detection bandwidth of 100kHz, 
we obtain a transmission record that averages over the 
effects of axial motion. 

A series of quantitative comparisons demonstrates that 
the remaining transmission signal faithfully reflects ra- 
dial motion with very little contamination from the axial 
averaging^. Simulations indicate that, for the parame- 
ter regime employed in Ref.Q, a trapped atom is typi- 
cally confined within ~ 50nm of a single standing-wave 
antinode until it heats very quickly and escapes the trap 
altogether. Events involving large-amplitude oscillations 
or 'skipping' across wells and retrapping in another antin- 
ode are very rare occurrences in this regime; since axial 
motion typically has such small amplitudes, we expect 
the average over it to have a negligible effect on measured 
cavity transmission. Comparison of observed maximum 
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FIG. 2: (a)First few levels of the ladder of energy eigenstates 
for the atom-cavity system. Uncoupled and maximally cou- 
pled eigenvalues are illustrated, along with their smooth de- 
pendence on the atomic position (and hence the atom-cavity 
coupling). Strength of driving light determines the system's 
population distribution across the first few levels of the lad- 
der, setting the shape and depth of the effective potential, (b) 
Effective potential for the atom, cavity, and probe detuning 
parameters of the simulations, and drive strength rihi corre- 
sponding to 0.3 photons in the empty cavity on resonance. 



transmission levels with the expected theoretical maxi- 
mum confirms this notion, giving an estimate of < 75nm 
for typical axial excursions. 

As further confirmation that transmission accurately 
reflects radial motion, we consider the observed period of 
transmission oscillations as a function of their amplitude. 
From our knowledge of the anharmonic (approximately 
Gaussian) trapping potential in the radial direction, it 
is straightforward to calculate the expected relationship 
between amplitude and period for these oscillations. As 
evident from Figure EHb), actual data closely follow the 
theoretical curve, indicating that the transmission record 
can indeed be interpreted as a record of radial motion 
(T(t) — > p(t)). Though the atom-cavity coupling is cylin- 
drically symmetric and thus provides no explicit informa- 
tion about 6(t), knowledge of pit) and the trapping po- 
tential allow us to reconstruct an estimate of an atom's 
angular momentum L(t) and thus of a two-dimensional 
trajectory in the (p, 9) plane. Such a method can be 
applied with success in a parameter regime where the 
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FIG. 3: (a) Cavity transmission record for an atom trapped 
with 0.3 photons in the empty cavity. Motion of the trapped 
atom within the cavity mode can be tracked via the oscilla- 
tions in cavity transmission. Each oscillation has amplitude A 
and period P as indicated on the trace, (b) Oscillation period 
P vs. amplitude A of transmission oscillations as in (a) . The 
solid curve is calculated for motion in the known anharmonic 
potential of Figure [^Jb), with no free parameters in the fit. 



atomic motion is largely conservative and the angular 
momentum varies slowly on the timescale of a single ra- 
dial orbit 0,0] ■ 

Three basic ambiguities will be clear from this algo- 
rithm for trajectory reconstruction: 1) the sign of the 
angular momentum is unknown, so the trajectory has 
arbitrary handedness. 2) the initial angle 8o is arbitrary, 
so the resulting trajectory can be rotated freely as a unit. 
3) The trajectory is constructed in two dimensions, with 
the axial motion confined within a single antinode, but no 
information is available about which antinode the atom 
occupies during the trajectory. These ambiguities, while 
noted here for clarity, arise in aspects of the motion not 
used in the feedback scheme treated below. 

Two-dimensional trajectory reconstructions dramati- 
cally illustrate the cavity-enhanced sensing power for 
atomic motion. However, the initial goal of our feed- 



back algorithms will be to control p(i); for this purpose 
it is sensible to ignore 8(t) and apply all available signal 
to noise to the task of estimating p{t) and p(t) in real 
time. The goal of such a program is then to use this in- 
formation to drive p(t) to a constant value, or in other 
words to circularize an orbit in the (p, 9) plane while not 
necessarily driving it to the cavity axis (p = 0). The 
latter task, which requires an explicit method of break- 
ing cylindrical symmetry for position sensing and for the 
effective potential, can be considered as a later extension. 



III. THE ATOM AND CAVITY AS A CONTROL 
SYSTEM: BASIC FEEDBACK STRATEGY 

As a guide in the identification of plausible feedback 
strategies and their limitations, it is useful to restate the 
problem somewhat in the language of control systems. To 
this end, we begin by setting aside the issue of axial mo- 
tion and treating the atom as a particle in a cylindrically 
symmetric, approximately Gaussian two-dimensional po- 
tential whose depth is controlled by the input light inten- 
sity: 



U fa -U e 



2 i 2 



(4) 



Note that the potential waist w is not simply equated 
with the previously introduced cavity field waist wq or 
with the mode intensity waist wq/V2, but rather is 
set by the self-consistent interaction of atom and light 
field in the strong coupling regime. Whereas the cavity 
mode profile is exactly Gaussian, U is only approximately 
Gaussian and has an exact form that is nonanalytical as 
determined by steady-state solutions to the master equa- 
tion for an atom at each value of p. The potential depth 
Uq depends on the intensity of the optical field used to 
drive the cavity mode. The potential waist w is in fact 
a (slowly varying) function of the drive strength as well 
0. 

The trapped atom is also subject to friction and to dy- 
namical noise (momentum diffusion), both arising from 
decays and re-excitations of the system on timcscalcs 
faster than the motion. In the regime of Ref. 0, the 
contribution of friction is small compared to the momen- 
tum diffusion terms in the equation of motion. 

Because the two-dimensional potential is symmetric, 
the atom's angular momentum L is constant, or rather, 
varies only due to dynamical noise. We can thus write a 
one-dimensional effective potential in the p dimension, 



V ef f = -U e- p2 ' 1 



L 2 /2mp 2 



(5) 



and thus an equation of motion (for an atom of mass m) 
2pU %-» 2 /™ 2 + -S, (6) 



P : 



m 2 p 3 



which we notationally simplify to the form 



-pEe~ 



(7) 
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where p is dimensionless (p — p/w), E — 2Uo/mw 2 is the 
input we control by varying the driving field strength, 
and L = L 2 /m 2 w 4 is constant except for the influence of 
friction and dynamical noise. 

The measurement of light transmitted through the cav- 
ity, T(t), is equivalent to a (noisy) measurement of pit). 
The noise of this measurement is largely fundamental 
quantum noise (shot noise) of detection. The mapping 
between T and p, derived again from steady-state solu- 
tions of the master equation for the coupled atom-cavity 
system, is not linear and furthermore depends on the 
value of the driving strength E. The initial objective is 
to circularize the two-dimensional orbit - in other words, 
to make p constant or to hold p = by varying the con- 
trol input E. 

The simplified system can be described by a block dia- 
gram as shown in Figure 0| The system exhibits myriad 
nonlinearities; for example, (T —* p) is nonlinear and 
depends on E, the dynamical noise depends on E, and 
the equation of motion for p is itself nonlinear. Nonethe- 
less, while this statement of the problem does not suggest 
provably optimal feedback strategies, it does motivate 
some conceptually simple algorithms based on switching 
between discrete values of driving strength E. Switching 
strategies of this sort are often invoked for the sake of ro- 
bustness, a major consideration in this scenario; robust- 
ness to dynamical and measurement noise is certainly 
important, but perhaps even more relevant is robustness 
to small uncertainties in system parameters (e.g., driving 
strengths or detunings). Switching or "bang-bang" type 
algorithms [2^] have the additional virtue of admitting 
easy implementation in simulations and in experimental 
design. 
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FIG. 4: Block diagram for the atomic position feedback loop, 
illustrating sources of noise and system nonlinearities. 

The feedback algorithms we investigate in this paper 
all share the same basic strategy of switching the driving 
field intensity between two discrete levels. This corre- 
sponds to switching between two potential depths (and, 



incidentally, two different sets of friction and momentum 
diffusion coefficients as well). The simple objective is to 
time this switching relative to the atomic motion so that 
an atom sees a steep potential when climbing out of the 
trap (p > 0) and a shallow potential when falling back 
towards the trap center (p < 0), as illustrated in Figure 
[5] The feedback algorithm, then, is based on switch- 
ing the potential at turning points of p, i.e., each time 
p crosses zero. Implemented effectively, this approach 
promises significant dynamical cooling of the radial mo- 
tion (p — ► constant or p — ► 0) in just a few oscillations. 
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FIG. 5: General feedback strategy for atomic radial coordi- 
nate. 



The initial detection and trapping of an atom are ac- 
complished as in Ref . 0, Q ; a weak probe at driving level 
exlo is used to detect the atom's arrival in the cavity with 
minimal effect on the motion, and an increase in trans- 
mission of this beam triggers a switch to driving level 
hi to populate strong-field-seeking states and trap the 
atom. Feedback is then implemented by switching the 
trapping potential between the hi level and an interme- 
diate lo setting, with switching times based on real-time 
information about the motion of the single atom. 

The simplest algorithm would be to switch back and 
forth between hi and lo potentials at the turning points 
of p(t), which are the zero-crossings of p(t). That is, trap 
initially in hi, switch to lo when p crosses zero from above 
(i.e., when p begins to decrease), switch back to hi when 
p crosses zero from below (p reaching its minimum and 
increasing), and so on until the atom escapes. However, 
this strategy calls for a theoretically infinite sequence of 
switching events, while it is desirable to instead achieve 
a steady state in some long-time limit. The presence 
of dynamical noise implies that the exact steady state 
of p — > is in any case unreachable, so we replace it 
with a goal of confining p to some range [— lim, +lim]. 
Thus the feedback strategy is modified to include slight 
hysteresis: lo—>hi when p — > -\-lim from below, hi—>lo 
when p — r —lim from above. With this modification, 
switching stops once p is confined within the acceptable 
range. Furthermore, we prefer a steady state with hi 
potential for reasons of the deeper confinement; to bias 
the system towards this final state, we use asymmetric 
hysteresis limits: lo^hi when p — > +lim from below, 
hi^lo when p —> —(lim + S) from above. 
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IV. SIMULATIONS OF FEEDBACK 
ALGORITHMS IN OPERATION 

For our simulations, we first choose driving strengths 
n-exio = 0.05 photons in the empty cavity, tim = 0.3 
photons in the empty cavity, and ri; = 0.15 empty-cavity 
photons. These driving strengths are high enough so that 
an atom of typical kinetic energy can be trapped by both 
lo and hi drives, yet low enough to ensure the increase in 
momentum diffusion between lo and hi does not outstrip 
the increase in potential depth. 

Our simulations of the atom-cavity dynamics are based 
on the treatment described in detail in Ref. 0, HH; 
the treatment is fully quantized for the atomic internal 
state and the cavity light field, but considers the atomic 
center-of-mass motion quasiclassically. This approxima- 
tion is suitable for the current experimental situation, 
with more manifestly quantized motion to be accessed by 
better cooling and/or detection of the atom's axial mo- 
tion. The non-conservative terms in the system, in the 
form of friction and momentum diffusion, are included in 
the simulation; the resulting "heterodyne transmission" 
trace is a perfect record of |(a)| 2 , on which measurement 
bandwidths and shot noise must be imposed separately. 
Shot noise is modelled as Gaussian white noise with an 
amplitude that depends on the size of the (noiseless) 
transmission signal. 

In the presence of sensor noise, we require estimators 
for the quantities p and p. Because one parameter (L) in 
the equation of motion is unknown and in fact slowly 
varying, we have chosen not to implement estimators 
based on a Kalman-filtcr approach. More sophisticated 
treatments include Kalman-type approaches to simulta- 
neously estimate p, p, and L, but these have not been 
explored in full detail. Meanwhile, we choose to estimate 
p(t) and p(t) directly from the measurement record, with 
no explicit reference to the equation of motion for the 
system. 

The noisy transmission signal T is sampled at 1MHz 
as in the experiment of Ref. |lj . To estimate p, we first 
perform an RC low-pass filter on T at 100kHz. This step 
is an infinite impulse response (IIR) filter which intro- 
duces only a small delay in the estimator. This filtered 
transmission signal is then put through a lookup table 
with linear interpolation to obtain p est (t). 

The resulting p es t(t) tracks the actual p closely but 
still with significant noise. Obtaining a time derivative 
without excessive noise thus requires some care. A va- 
riety of methods are mentioned in Ref. 30], in which 
the authors are concerned with estimating the sign of a 
time derivative in order to feed back to a system - es- 
sentially the same problem we encounter. We employ a 
simple finite impulse response (FIR) filter that takes the 
slope of a linear least squares fit to p es t (t) over a window 
of fixed size. A detailed implementation of this filter is 
found in Ref. 0. The resulting p est (t) is a good esti- 
mator for p at the middle of the window, so the delay 
induced is approximately half the window size. We find 



via numerous simulations that a window size of 30 to 40 
ps gives a signal p es t(t) which is quiet enough for use in 
our control. Thus reliable estimation in the presence of 
noise introduces a delay of approximately 15 to 20 ps in 
the feedback loop. This time delay can be compared to a 
typical atomic orbital period of r r ~ 100/xs, correspond- 
ing to a period of ~ 50ps for p. Feeding back effectively 
in the presence of such large delays requires a certain 
amount of adjustment to the naive cooling algorithm, as 
discussed in detail below. 



A. Actual Dynamics but No Measurement Noise 

Before treating the case of actual experimental noise, 
we explore the performance of our feedback strategy in 
simulations with noiseless measurement and thus perfect, 
zero-delay sensing of (p, p) . Figure shows an exam- 
ple trajectory using this asymmetric-hysteresis switching 
strategy. The values of cavity transmission T and atomic 
position variables are sampled every 1/Lts, but the dynam- 
ical timestep is 3,000 times finer than this "information" 
timestep. Note that axial motion (the x direction our 
strategy neglects) is included in the simulation, and when 
its amplitude is large it gives rise to the very fast varia- 
tions seen in T(t). However, since the period of x motion 
is similar to the information timestep used, note that 
these signals are undersampled in the record. 
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FIG. 6: Simulated feedback example with perfect sensing of 
p(t). The graph on the left shows idealized cavity transmis- 
sion T(t) (purple), radial position p(t) (blue), and dp/dt (red) 
for a simulated atom trajectory. In the right panel the same 
atomic trajectory is shown in the (x, y) (or (p, 6)) plane. This 
trajectory begins as a nearly vertical transit through the cav- 
ity, but is circularized by the feedback employed. In the left 
panel, the circularization is evident as dp/dt is damped to- 
wards zero and p and T become nearly constant at the end of 
the trajectory. 

A 10-ps box filter is applied to p in order to remove 
some oscillations caused by x motion and also partially 
to anticipate some effects of noisy detection and de- 
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lay. In setting the conditions for potential-switching, 
we employ the asymmetric hysteresis described above 
with (lim, 5)=(0.05,0.03)/xm//is, so the potential depth 
is switched at p — +0.05, — 0.08^m/^s. Switching events, 
since they correspond to turning the light level up and 
down, can be seen as sharp edges in the transmitted light 
T. As the example illustrates, the control strategy suc- 
cessfully circularizes atomic trajectories within a few or- 
bital cycles. This can be seen from p(t) as well as from 
the trajectory shown in the (p, 9) plane. The hysteresis 
limits are chosen so that variations in p due only to dy- 
namical noise tend not to trigger any switching of the 
drive. This is illustrated by the continued high control 
level throughout the time t = 500 — 650/is, while p is wan- 
dering diffusively rather than oscillating with regularity. 

The overall trap lifetime is dominated, as in this exam- 
ple, by heating in the x direction. (In the example shown, 
note the fast, large-amplitude variation in transmission 
just before the atom escapes; this is a signature of rapid 
axial heating.) Thus our feedback strategy has little im- 
pact (at the level of 10%) on average trapping lifetimes. 
Circularizing the orbit helps decrease axial heating since 
the potential depth no longer wanders as p varies; how- 
ever, the feedback is accomplished by sharp switching 
events which occur at arbitrary times relative to the oscil- 
lations in the x direction. The overall impact on lifetimes 
is therefore small in the simulations we have performed. 
Since the feedback algorithm is aimed at reducing mo- 
tion in the p direction, its success is best measured by its 
performance at that task specifically. Lifetime effects can 
become apparent only if the axial motion is suppressed by 
some other means; that case is treated briefly in Section 
E below. 



B. Adding Measurement Noise Adds Delays 

The addition of measurement noise and consequent es- 
timation of p introduce significant loop delays, as de- 
scribed above. Since p es t can be almost half a cycle be- 
hind the actual p, we expect naive switching to be well 
out of phase with the atomic motion and thus relatively 
ineffective as a cooling mechanism. Figure[7|shows an ex- 
ample in which the feedback strategy is identical to that 
of Figure El but applied to p es t rather than to p itself. 
The resulting time delay seriously compromises perfor- 
mance, as shown. In the figure, note that the switching 
events, recognizable as sharp edges in transmission, do 
not line up with turning points of p. As a result, p is not 
damped and the trajectory remains elliptical. 



C. Account for Delays by Waiting a Cycle 

Since it seems clear we cannot simply close the loop 
with the delays necessitated by measurement noise, we 
choose to address the problem by adding even more delay 
- that is, by detecting a switching condition (p es t crossing 




100 200 300 400 500 600 700 
Time (|is) 



Trajectory in (p,6) plane 




Horizontal Position (jim) 



FIG. 7: Measurement noise leads to filtering and unaccept- 
able loop delay in this example. In addition to the quantities 
plotted in Figure|S| the left panel now displays the experimen- 
tally realistic transmission T(t) with noise (yellow), and the 
radial motion estimators p e st(t) (green) and (dp/dt) e3t (light 
blue) arising from the noisy transmission signal. Feedback 
switching is triggered by the value of (dp/dt) ea t, but since 
the estimator lags well behind the true dp/dt damping is not 
achieved, and indeed the two-dimensional trajectory shows 
little change in shape over time. 



a hysteresis limit) and then waiting to switch the poten- 
tial at a time which should catch the next oscillation of 
p. A first attempt in this direction would be to assume 
a fixed period for oscillations of p. In this case the addi- 
tional wait before switching is given by this fixed period 
minus the estimator delay for p es t- Since each switching 
time is now set by the detected signal from the previous 
cycle, the first switching opportunity (first minimum of 
T and maximum of p) will be missed in this strategy. 
Rather than miss this cooling cycle, we impose a single 
switching event a fixed time after the initial trap turn- 
on. Thus the potential switches exlo^hi on the initial 
trigger, hi—>lo a fixed time later, and lo<-^hi thereafter 
based on the last zero-crossing time of p es t ■ 

However, the actual dynamical period varies by eas- 
ily a factor of two over the course of an atom's trapping 
lifetime due to changing amplitude of oscillation in the 
anharmonic potential, as seen in Figure 3b. Thus the 
fixed-period assumption is a poor one. A better strat- 
egy is to record the length of each period in p est and 
assume each cycle will be the same length as the previ- 
ous recorded one. Thus the "waiting time" estimate will 
adjust itself as the dynamical period changes, though it 
will in general be one cycle behind. This strategy is em- 
ployed for the trajectory shown in Figure |S] The initial 
switch occurs 45/iS after trap turn-on, the least-squares 
window is 40/is, and the "wait time" between subsequent 
p es t limit-crossings and the resultant potential switches is 
given by the previous period minus 20/is. This switching 
strategy, with deliberate delay based on an active mea- 
surement of the p oscillation time, appears to be a viable 
means of performing control in the presence of sensor 
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noise and its associated loop delay. 




Trajectory in (p,9) plane 




Horizontal Position (pm) 



FIG. 8: In this example delay is dealt with by tracking p 
turning points and using this information from each cycle to 
switch the potential at the predicted next turning point. This 
strategy give moderate success at damping dp/dt and circu- 
larizing the trajectory. 
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D. Comparisons with Open Loop Strategies 

To evaluate the effects of feedback more quantitatively, 
we introduce a figure of merit for the damping of radial 
oscillations in an atomic trajectory. Since the goal of 
the control strategy is to confine p near zero, the per- 
formance can be measured by comparing the variance of 
p over intervals of equal duration near the beginning of 
the trajectory and after feedback has been operating for 
some time. We choose a time window of duration 200^s 
as long enough to encompass well over one cycle of the ra- 
dial motion. The comparison is taken between two such 
windows separated from one another by 200/is; this delay 
is selected as long enough for several cycles of feedback 
action, yet short enough so that the statistic exists for a 
large fraction of trapped atom events. Finally, we base 
our performance measure on the experimentally accessi- 
ble p es t rather than on p itself. Thus our figure of merit 
for feedback performance is given by 



M = 



(Pest(t)) 
(f>est(t)) 



(8) 



where times are measured from the initial trapping 
(exlo^hi) switch. Large values of the quantity M cor- 
respond to well-damped radial motion, p(t) — > constant, 
though orbits may still be circular at any radius p > 0. 
(Damping in the sense of actual energy removal is dis- 
cussed explicitly in Section E.) Small (~ 50 — 100/is) 
changes in delay time or window size have been investi- 
gated and do not appreciably change the nature of the 
results for M. 

FigureElshows histograms of M for several data sets in 
which different switching protocols, detailed in Table [I] 




open loop (01,02) 



(e) 



open loop (03,04) 
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FIG. 9: Figure of merit for feedback in closed- and open- 
loop cases (full dynamics). M >1 indicates damping of radial 
motion. 



have been employed. Each data set is generated by sim- 
ulating a fixed number of individual atom drops from the 
known distribution of initial conditions. Only some frac- 
tion of trajectories result in a triggering/trapping event, 
and of these only a fraction of atoms have dwell times 
long enough to compute a value for M. Thus, for exam- 
ple, set C2 was generated from 5,000 trajectories, yield- 
ing 1,335 trigger events and 147 trajectories for which M 
could be obtained (i.e., with dwell times at least 615/is). 

While Table [I] gives the specifics of each data set rep- 
resented in Figure El the essential comparison is between 
closed loop - i.e., active feedback - algorithms and open 
loop counterparts which simply switch potentials in a 
predetermined sequence independent of real-time posi- 
tion information for the individual atom. The closed loop 
algorithm is that of Figure|Hl in which measurement noise 
and loop delays are dealt with by waiting nearly one cycle 
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Data set 


switching protocol 


# atom drops 


# trapped 


CI 


Closed loop 


2000 


534 




Switch once 








45 ps after 








initial trigger. 








Thereafter switch 








(previous cycle length) 








- 20^is after 








(dp/dt)est 








crosses limits. 






C2 


Same as above (CI) 


5000 


1335 


01 


Open loop 


2000 


561 




Switch every 45 ps 








after initial trigger. 






02 


Same as above (Ol) 


5000 


1319 


03 


Open loop 


2000 


552 




Switch every 35/is 








after initial trigger. 






04 


Same as above (03) 


5000 


1322 



TABLE I: Exact conditions used for data of the full- 
simulation histograms. 



ent regime from the one in which our simulations operate. 
In particular, while both scenarios involve a measurement 
bandwidth which averages over axial motion, in the case 
of Ref. 01 the amplitude of the axial motion is in fact 
quite large. One key observation in that experimental 
setting was in fact the change in measured transmission 
as atoms either oscillated within a single standing-wave 
antinode or "flew" across multiple antinodes. Without 
the measurement bandwidth to observe axial oscillations 
directly, one is hard pressed to separate axial modulation 
from radial motion in the manner presented here. While 
both feedback algorithms involve a similar switching of 
cavity driving strength, the feedback modeled here is di- 
rected at cooling a particular component of the atom's 
motion r(t). Figures l§l and ITUI illustrate the damping of 
dp/dt. By contrast, the lifetime enhancement of Ref. 17] 
occurs because the authors are able to discern when the 
atom is trapped at high g(r) and selectively turn down 
the otherwise large diffusive heating by lowering the trap 
intensity at those times. 



to apply the knowledge of motion gained during the pre- 
vious oscillation. The open loop algorithms, in contrast, 
simply switch the potential between hi and lo at fixed 
intervals following the initial trapping event; the fixed 
interval is chosen to coincide with a reasonable average 
value for an atomic oscillation period. Over a long atomic 
trajectory, the atomic trajectory clearly evolves out of 
phase with any one open loop switching algorithm; since 
this evolution is different for each atomic trajectory, we 
compare closed-loop with open-loop strategies as a means 
of evaluating the importance of real-time measurement 
and feedback in our algorithm. 

Closed loop, active feedback clearly damps radial os- 
cillations more effectively than its open loop counterpart. 
The M distributions in Figure H0[a,c) have larger mean 
than those in Figure 0Jb,d,e), with the closed-loop his- 
tograms showing many trajectories pushed out to higher 
values of M by the active feedback. These results indi- 
cate that real-time measurements of p(t) can indeed be 
applied to facilitate cooling of a single atom's motion in 
that dimension. 

Further refinements should improve the performance of 
the algorithm. For instance, the cycle-length predictor 
could be changed to allow asymmetries between p > 
and p < half-cycles. Additionally, the filters them- 
selves could be adjusted or replaced with better estima- 
tors which incorporate information about angular mo- 
mentum L. 

The simulations reported here could in principle be 
employed, with minor modifications, to address the ex- 
perimental regime of the atom-cavity feedback performed 
in Ref. 0| . We have not attempted a quantitative com- 
parison since results may be sensitive to details of the 
experimental implementation. However, we note that 
because of the values of (<?, k) for the optical cavity em- 
ployed, Ref. [13 was carried out in a qualitatively differ- 



E. Performance with Axial Motion Suppressed 

In a final set of simulations, we investigate the perfor- 
mance of our radial cooling algorithm in a setting where 
the axial atomic motion is independently suppressed. 
With no (or little) axial motion, axial heating no longer 
limits trapping times and the effects of radial feedback 
can be seen more clearly. To achieve this in simulations, 
we impose an ad hoc elimination of the axial dimension; 
however, this case could be relevant to several future ex- 
perimental scenarios. For example, trapping and sens- 
ing mechanisms, both currently accomplished with the 
same probe beam, could be separated to allow a trap- 
ping field with a low scattering rate and much-reduced 
axial diffusion. Alternatively, the separation of axial and 
radial timescales could be exploited, either to apply axial 
cooling 32] between cycles of radial feedback or to simply 
avoid extra axial heating by ramping the potentials up 
and down at a rate that appears adiabatic to the axial 
motion while still impulsive in the radial dimension. 

Prospects for implementation aside, simulations with 
no axial motion demonstrate some aspects of the radial 
feedback protocol that are otherwise less transparent. We 
explore this regime with a set of simulations that differ 
in three ways from those presented above. First, the 
axial dimension is eliminated entirely and the atom is 
artificially constrained to remain at rest at an antinode 
of the standing-wave cavity field. Second, since axial 
heating is no longer an issue, we employ somewhat deeper 
trapping potentials than in the simulations above. The 
weak probe level is still n ex i a = 0.05 photons in the empty 
cavity, but now we turn on the trap initially at a level 
n ex hi — 0.6 photons in the empty cavity, and we feed back 
by switching between this and the weaker level nhi = 
0.3 photons. The effective potentials thus generated are 
~ 50% deeper than in the previous simulations using 
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nhi and rii . Finally, without axial motion we employ a 
coarser computational timestep of (l/30)/xs. 

Figure EJa) shows the feedback figure of merit for the 
cases of closed-loop feedback, constant trapping at nhi 
or n ex hi, and open-loop switching. The open- loop proto- 
col in this case is to trap initially with n ex hi and switch 
between driving levels n ex hi and nhi every AOps during 
the transit duration. Taking advantage of longer over- 
all atom dwell times, we can now consider time windows 
separated by a greater delay, so the quantity displayed 
here is 



(a) 

Simulation feedback performance for measurement noise but no axial motion 



M' 



{Pest{t)) 
' 1015 us ->1215 fJ.s (Pest(t)) 



(9) 



rather than the original M of Figure The result is 
qualitatively similar to that of the full simulation, with 
the closed loop strategy performing significantly better 
than its open loop counterpart. In this case, the mean 
value of M' for closed-loop feedback is ~ 2.5 times greater 
than for open-loop switching or for no switching. 

The data of Figuresl^land lTUT a) indicate that our feed- 
back algorithm acts to drive p to zero, i.e., to circularize 
atomic orbits at a constant value of p. If angular mo- 
mentum L is not correspondingly increased, this effect 
implies a damping of radial energy. However, we may in- 
vestigate more directly whether this algorithm actually 
removes total energy from the radial motion. With ax- 
ial motion eliminated, we can now explore this issue by 
asking how the feedback algorithm affects trapping life- 
times. Figure lTUT b) shows atom dwell times for the same 
three cases of closed-loop feedback, constant drive levels 
of fiexM or nhi, and open-loop switching. The increase in 
lifetime for closed-loop feedback is immediately apparent. 
Indeed, the closed-loop results agree well with an expo- 
nential lifetime of 8.9ms, as contrasted with 2.6ms for 
trapping at n e xhi alone, 1.9ms for trapping at nhi, and 
1.1ms for open-loop switching between the drive levels. 

In the closed-loop case the trap potential is varied dur- 
ing the transit but never exceeds the depth associated 
with driving level n ex hi- Nevertheless, trap times ex- 
ceed those for constant driving at n ex hi, demonstrating 
that active feedback as applied to the radial dimension 
does act to remove radial energy, in addition to simply 
pinning p to a constant value. The same point can be 
illustrated by considering the change in total energy be- 
tween the beginning (15-215 ps after trigger) and end 
(1015-1215 ps after trigger) of an atom's dwell time in 
the cavity; the active feedback strategy produces a mod- 
est (~10%) net energy removal not seen under either the 
simple trap or the open-loop switching protocol. Since 
the feedback algorithm performs better than both open- 
loop and fixed-trap strategies when measured by dwell 
time, radial damping M', or total energy removal, we 
characterize its effect as in fact actively cooling a com- 
ponent (p{t)) of the atomic motion. 

Note that all lifetime in the two-dimensional simula- 
tion are enhanced relative to the full three-dimensional 
case, in which both experiments and simulations have 
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FIG. 10: Feedback performance and lifetime enhancement in 
simulations with no axial motion. M' >1 indicates damping 
of radial motion. 



mean trapping times of only ~ 400/is. Lack of lifetime 
enhancement from radial feedback in the full simulation 
is some indication of the very weak mixing between axial 
and radial motion, so that cooling in one of these dimen- 
sions does little to control temperatures in the other. 



V. OUTLOOK FOR EXPERIMENTAL 
IMPLEMENTATION 



The feedback simulations discussed in this paper have 
been conducted with very close reference to the condi- 
tions realized in the experiment of Ref. p] , in particular 
for cavity properties, trapping statistics, and signal-to- 
noise in the balanced heterodyne detection. The cur- 
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rent experimental effort, aimed at realizing the feedback 
strategies described here, employs very similar condi- 
tions while incorporating some improvements as in Ref. 
LH ■ Notable changes from Ref. Q include a slightly 
shorter cavity, improved vacuum pressure in the cavity 
region (3xl0 -10 torr) enabled by a differentially pumped 
double chamber, and cavity length stabilization via an 
error signal generated by an independent laser one free 
spectral range away from the cavity QED light. Imple- 
mentation of the feedback strategy described above is 
clearly outside the regime of analog electronics, so an 
additional modification is the use of digital processing 
and FPGA technology |33|. With these tools, experimen- 
tal data similar to the simulation results presented above 
seem well within reach. 

It seems reasonable to ask how much experimental data 
should be necessary to exhibit a distinction between ac- 
tive feedback and open loop schemes. From the simula- 
tions of Figure (Table |]J, we see that with data sets of 
about 500 trapped atoms the differences between open- 
and closed-loop schemes already begin to become appar- 
ent, and these differences are well demarcated with two 
or three times that much data. With fairly conservative 
estimates of one trigger per MOT drop and one MOT 
drop every 5 seconds, this means significant effects could 
well be seen with just one to two hours of experimental 
data at each setting. Much more data collection is ex- 
perimentally realistic, allowing exploration of a wealth of 
additional questions. 

With the atom-cavity system's capacity to give real- 
time information on a trapped atom's position, active 
feedback might seem to be an attractive enabling tech- 
nology for experiments that require a stationary atom 
at fixed g(r) in the cavity field. The present work was 
undertaken in part to explore this option through simu- 
lation of realistic experiments. While we conclude that 
feedback of a uniquely real-time nature can measurably 
alter properties of the atomic motion, we also find that 
measurement bandwidths and signal to noise strongly 
constrain the precision of feedback performance under 
current experimental conditions. Feedback experiments 
address important topics in quantum measurement and 
control; meanwhile, trapping strategies using auxiliary 
fields offer a more direct route towards a stationary atom 
in a cavity. 



VI. CURRENT LIMITS AND FUTURE 
DIRECTIONS 

The feedback algorithm developed above for the atomic 
radial position p is subject to basic limits arising from 
dynamical and measurement noise in our system. These 
limits can be expressed as lower bounds on the (one- 
dimensional) temperature for p. Because the strategy is 
based on discrete switching, with feedback delayed and 
timed based on the previous switching-cycle length, the 
control is always based on information gathered over the 



previous motional cycle. Thus dynamical noise over an 
atomic motional timescale will set a lower limit on T p . 
Referring to , we find that momentum diffusion (due to 
spontaneous emission) gives a typical energy increase per 
radial oscillation time of AE Tr ~ 0.02([/ ) ~ ksi^OpK). 
Furthermore, measurement noise places a limit on the 
detectable amplitude of p variations. This amplitude de- 
pends strongly on the absolute value of p due to the non- 
linearity of the T — > p mapping. However, using the mea- 
sured sensitivity from the atom-cavity microscope, we es- 
timate that over a motional cycle we can resolve p oscilla- 
tions of amplitude (20nm/^jHz)^l/2iTT r « 0.77/im. On 
the side of the cavity mode, where the effective potential 
is steepest, this corresponds to T p sa 150pK . While this 
limit corresponds closely to the simulations of the pre- 
vious section, where axial motion is suppressed, the full 
simulation never reaches this limit because axial heating 
cuts off atomic lifetimes too quickly. Thus improvements 
to address axial heating are of great interest for seeing 
the full effect even of radial cooling. 

Beyond the experimental and algorithmic variations al- 
ready discussed, a number of broader questions are raised 
by the use of active feedback to dynamically cool a com- 
ponent of motion for a single atom. One question deals 
with the ultimate limits of such a cooling mechanism. 
Within the current experimental setting, limits to radial 
cooling arise from atomic lifetimes (dominated by axial 
motion) , but are also constrained by the dynamical noise 
and by the shot noise of detection. Some lifetime and 
dynamical noise issues could be addressed by separating 
trapping and sensing, for example by using a far off reso- 
nance trap in conjunction with a sub-photon level cavity 
QED probejj]. The remaining issues would then center 
on signal-to-noise for the atomic position measurement, 
as well as on limits imposed by backaction of the mea- 
surement itself as it approaches the standard quantum 
limit [iHI^ . These limits must be considered not only in 
the context of near- resonant probing, as treated here, but 
also in the case of a far detuned probe for which atomic 
position information is extracted from probe phase shifts, 
as first measured in [20| . 

Extension of active feedback beyond the p dimension 
raises related questions. The question here is one of using 
various techniques - for example, a symmetry-breaking 
potential as could be provided by a higher-order trans- 
verse mode of the cavity or frequency-domain filtering of 
the signal - to estimate and control a three-dimensional 
vector using the time record of a single quantity, the 
transmitted light field. Undoubtedly, different driving 
parameters, detection methods, and data processing will 
be appropriate depending on the relative importance 
placed on information about each dimension of the mo- 
tion. These questions address in various ways some basic 
issues of optimal state estimation and control for single 
quantum systems, and this experimental system promises 
to be a rich one for exploring such issues in further detail. 
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